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Abstract — We study amplifled-and-forward (AF)-based two- 
way relaying (TWR) with multiple source pairs, which are 
exchanging information through the relay. Each source has 
single antenna and the relay has multi-antenna. The optimal 
beamforming matrix structure that achieves maximum signal- 
to-interference-plus-noise ratio (SINR) for TWR with multiple 
source pairs is derived. We then present two new non-zero- 
forcing based beamforming schemes for TWR, which take 
into consideration the tradeoff between preserving the desired 
signals and suppressing inter-pair interference between different 
source pairs. Joint grouping and beamforming scheme is pro- 
posed to achieve a better signal-to-interference-plus-noise ratio 
(SINR) when the total number of source pairs is large and the 
signal-to-noise ratio (SNR) at the relay is low. 

Index Terms — Analogue network coding (ANC), two-way 
relaying (TWR), multiple source pairs, information exchange, 
analogue relaying, optimal beamforming. 



I. Introduction 

By applying the physical-layer network coding [l] or 
analogue network coding (ANC) ^ in two-way relaying 
(TWR), only two time slots are required for one complete 
information exchange using TWR. In the first time slot, both 
source nodes transmit simultaneously to the relay. In the 
second time slot, the relay broadcasts the common message 
which is obtained by combining the received messages. Since 
both source nodes know their own transmitted signals, each 
of their self-interference can be completely canceled prior to 
decoding. 

The TWR has been studied in lO to ifTOl for the case 
of single source pair. The beamformer design for AF-based 
multiple-input-multiple-output (MIMO) TWR is studied in 
Sllll, in which the receive and transmit beamforming are 
derived separately and then combined to form the relay 
beamformer Furthermore, some strategies to enhance the 
performance of TWR can be found in ifTTl to ifTsl and 



references therein. 

The non-ANC-based TWR for multiple source pairs is 
studied in lfT4l to ifTSl . Unlike the single source pair case, 
in multiple source pairs scenario, additional inter-pair inter- 
ference exists between different source pairs, which degrades 
the TWR performance. In lfT4l . a relay network with multiple 
source pairs and multiple relay nodes is studied, where all 
sources and relay stations have multiple antennas. The mul- 
tiuser TWR is proposed and studied in |15|, where multiple 
source pairs are communicating via multiple relays. In iflTl 
|,18J , the MIMO TWR where multiple wireless node pairs are 
communicating via a single decode-and-forward (DF) relay 
is studied. 

Due to the presence of inter-pair interference, previous 
beamforming solutions for TWR with single source pair is 
no longer useful and new solutions are required for the case 
of multiple source pairs. In almost all the above works lfT4l 
to 1 18 1, the inter-pair interference are canceled using zero- 
forcing (ZF)-based approach. 

In this paper, two new non-ZF-based beamforming 
schemes or beamformers are proposed for ANC-based TWR. 
Instead of completely canceling the inter-pair interference 
for all source pairs by ZF-based methods as was done in 
previous works, we propose joint grouping and beamforming 
scheme that divides a given large number of source pairs into 
smaller subgroups, and then apply the proposed beamformers 
to each subgroup. To the best of our knowledge, this approach 
has not been studied in previous works. Simulation results 
are presented to compare the performance of the proposed 
schemes. 



II. System Model 

We study wireless TWR with a multi-antenna relay and 
a total number of Kt single-antenna sources, where Kt 
is an even number Due to the fact that the inter-pair 
interference contains the desired signals of all the other 
sources, these desired signals are also suppressed by any 
suboptimal beamformer which causes significant loss in the 
SINR, especially when Kt is large. We propose to overcome 
that this shortcoming by first dividing a large number of Kt 
source pairs into N subgroups, each with a smaller number 
of K source pairs, where K is an even number. Then, by 
using time division approach, the relay performs non-ZF- 
based beamforming on each subgroup of users at one time, 
and take turn to serve all source pairs, to achieve a better 
throughput performance. Next, we consider a given subgroup 
of K sources, and derive TWR beamformers for these K 
sources. 

Without loss of generality, we assume the fc-th source 
node Sk, k — . . . ,K, is to exchange information with 
another source S*^, where k = 1,...,K, k ^ k. Each fc- 
th source have single antenna, whereas the relay station R 
is equipped with M{M > K — 1) antennas. Let the M x 1 
vectors h^: and denote, respectively, the channel response 
matrices from Sk to R, and that from Sf. to R. We assume 
that the elements of and follow the distribution of 
circularly symmetric complex Gaussian with zero mean and 
unity variance, which is denoted as CAf{0, 1). Let Sk{n) and 
Sj,{n) denote respectively, the transmitted symbols from Sk 
to Sf., and from Sf. to Sk- We assume the optimal Gaussian 
codebook is used at each Sk, and therefore Sfe(n)'s are 
independent random variable each is distributed as CJ\f{0, 1). 

Assume TDD is used and two time slots are needed 
for information exchange using analogue network coding 
(ANC)|[T|, ||2l. In the first time slot, all active source nodes 
transmit their signals simultaneously, the received baseband 
signal vector y_R(n) at R is given by, 

ynin) =hi^ si{n) + ... + hK-i^/Pjf_^SK-i{n) 

(1) 

where n is symbol index, pk is the transmit power at Sk, 
Zfi{n) is the received additive noise vector, and without loss 
of generality, it is assume that Zfl (n) follows the distribution 
of CAf{0, Icr^), Vn, where I is an identity matrix. Throughout 
this paper, we assume that the p^'s are given or fixed. 



At the relay, we consider AF relaying using linear beam- 
forming which is represented by a M x M matrix A. 
The transmit signal x/j(n) at R can be expressed in terms 
of its inputs y_R(ri) as X]i{n) = Ay ]i{n). We assume 
channel reciprocity for uplink and downlink transmission 
through the relay. In the second time slot, when Xfl(n) is 
transmitted from R, the channels from R to Sk become h^, 
k — 1, . . . , K. The total transmit power at R, denoted as pn, 
can be shown as, 

K 

PR = Y^ liAhfell V + Tr(AAt)a2, (2) 

k=l 

where Tr(X) denotes the trace of X. 

The ANC is adopted as follows. We assume that using 
training and estimation, h^Ah^ and h^Ahj, are perfectly 
known at 5*^, fc — 1,...,K prior to signal transmission. 
Each of the Sk can first cancel its self-interference and then 
coherently demodulate for Sf,. This yields 

=hfc Ah^7^^s^(n) + J2 hfcAhjVPjSjW 

j^kj^k 

+ Zk{n),k ^l,.. .,K, 

for k = 1, . . . ,K . We assume that the received noise Zk{n) 
is distributed as CA/'(0,ct^), and ZTj{n) are independent of 
Zk{n). Zk{n) — h^Azfl(n) + Zfc(n), and Zk{n) is distributed 
as C7V(0, (||hjA|p + l)cr^). At each Sk, coherent signal 
detection can then be used to recover s,f.{n) from yk{n). 
The signal-to-interference-plus-noise ratio (SINR) for the fc- 
th destination node, k ~ 1, . . . ,K, can be expressed as 



' E,^fc,,^fc|hrAh,Pp,-(n) + (||h^A||2 + l)a2- 
III. Proposed Schemes 

A. Optimal beamformer 

We define the uplink (UL) channel gain matrix H[/i = 
[hi, h2, . . . , hA'] and denote the singular value decompo- 
sition (SVD) of HuL as Hc/L = USV^, where the 
M X K matrix U, and the K x K matrix V, are with 
orthogonal column vectors, and S is singular value matrix 
with dimension K x K where [S]^.^ = cr^ for £ = 1, . . . ,K. 
We have the following new result. 

Proposition: The optimal beamforming matrix to achieve 
maximum SINR in (O has the following structure: 

A°P' = U*B°P'U^, (4) 



where B°p* is ci K x K matrix. 

Proof: The above new result can be proven by extending 
the proof of |[8l|[lOl, which considers the TWR for the 
case of single source pair who are exchanging information. 
For the case of multiple source pairs, each source is also 
subject to interference from the other source pairs. This 
so-called inter-pair interference term also depends only on 
B as E,^fe,,^fc|hrAh,f = |hi^U*BU^h,f , 

which also spans the total signal subspace of B. Therefore, 
we obtain B°p*. The beamforming matrix B can be solved 
as shown next. 

Let hfe = U^hfe, k — 1, . . . ,K, represent the effective 
channel from Sk to R, with A given in (|4]i. Similarly, let 
represent the effective channel from R to 5*^. The SINR 
formula in Q can be written in terms of B. Throughout this 
paper, the optimization problems are formulated in terms of 
the beamforming matrix B and the effective channels h^. 

The minimum power (MP) beamformer is derived by 
minimizing the total relay transmit power with respect to the 
relay beamforming matrix B (or equivalently. A), subject to 
SINR constraints 7fc > 7^, fc = 1, . . . , -ft'. 



interference minimization problem is formulated as. 
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^|lBhfc|lV+Tr(BB^)a2, 



fc=i 



(5) 



subject to > Jk^k = 1, . . . ,K. 



For given transmit powers pk, k = 1,...,K, we can 
use the same approach as in ifTOl . to develop an efficient 
algorithm based on the second-order cone programming 
(SOCP) lHU to solve the problem in ©. 

B. Suboptimal beamformer 

We derive the suboptimal minimum interference (MI) 
beamformer that does not require computations via optimiza- 
tion technique. The tradeoff between the desired signals and 
interference is taken into account by minimizing the sum 
of inter-pair interference plus AF noise at the output of the 
beamformer, subject to the constraints that the desired signal 
gain for each fc-th receiver is equal to a constant /3fc. The 
additive Gaussian noise Zk{n) at each fc-th source, which is 
not affected by the beamforming matrix B, and is neglected 
in the optimization. Again in this case, the joint design of 
both receive and transmit beamforming is considered. The 
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subject to Ih^Bh^pp^ 



/3fe,fc = l, 



where 
4.(B) 



hlBh,|2r 



|hi^B||V,fc = l, 



(6) 



(7) 



The above constraints in Q are introduced to preserve the 
desired signal components at each source, so as to minimize 
the inter-pair interference plus AF noise component that point 
into any undesired direction. In (|7]l, //c(B) represents the sum 
of inter-pair interference (from other source pairs) plus AF 
noise that is imposed on Sk- To solve for B, the UL and 
DL channel response of each fc-th source is assumed to be 
known at the relay. 

We assume that the p^'s are given, and write 

K,bp, and ||hjB||2 = ||G,b|p, where fc ^ fc, 
and b denotes a x 1 equivalent beamforming weights 
vector which is generated by the rule of (|8]l as. 
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(8) 



The problem (|6]l can be written in terms of b as. 
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(9) 



subject to |f^ b| = ^fc, fc = 1, . . . , K. 
We further define Rfe = Ej^fcj^fed*^.d 



G^Gfc, where both R^, and are K"^ x if ^ matrices. By 
writing ||Gfcb||2 = b^Nfcb, E,^k.,^k = b^Rfcb, 

and # = ^^^(Rfc+Nfecr^), the problem (|9]l can be written 

as, 

min b^$b, 

subject to C b = g. 

where C = [fj* f| fg f^] is the K'^ x K constraint matrix, 
and ffc has been defined in the paragraph after Eq.®. Each 
of the column of C imposes a constraint on the equivalent 
beamforming weight vector b. The K x 1 response vector 
g = [/3i /32 • • • I^k]'^ contains the K scalar constraints values 
/3fc,fc = 1,...,K, which are chosen to satisfy the SINR 
requirements of each source. 



The beamforming vector solution, denoted as bMi, that 
satisfies ( fTOl i is solved as ll20i . 

bMi = *"'C(C^*-iC)-ig. (11) 

The MI relay beamformer, denoted as Bmi, can then be 
obtained from bMi as Bmi = ^^^(abMi), where V^^ 
denotes the inverse operation of V defined in (O. It can be 
shown that the solution Bmi maximizes the SINR of each 
source, where the constant a is used to control the total 
relay power p/j in (|2]). The iterative steps to search for a 
are presented next. 

C. Joint grouping and beamforming scheme 

We propose the following joint grouping and beamforming 
scheme that divides a given large number of source pairs into 
smaller subgroups, and then apply the above beamformers 
to each subgroup. For simplicity, here we consider that the 
grouping is done arbitrarily. By doing so we can reduce 
the feedback of channel state information and beamforming 
calculation during each relaying. 

• Divide the total number of Kt source pairs into N 
subgroups, N = 1,2,..., and let K = Kt/N denote 
the number of source in each subgroup. 

• For each n-th subgroup, n = 1, . . . ,N, compute the MI 
beamforming matrix as follows: 

. Given a e [0, amax], P = [pi P2 • ■ • Pk]'^- 
. Initialize aiowcr = 0, auppcr = "max- 

• Compute the SVD of the UL channel matrix Hul = 
USV^. 

• Estimate the correlation matrix 

• Compute the constraint matrix C and the response 
vector g. 

« Compute the beamforming weights using bMi = 
*-iC(C'f^*-iC)-ig. 

• Obtain the MI relay beamforming matrix by arranging 
the beamforming weights of bMi into the matrix form, 
Bmi = V^^iahui)- 

. Repeat 

1. Set a <- 5 (aiowcr + auppcr)- 

2. Update a by the bisection method fi9]: If pn is 
less than a given power constraint, set aiowcr ^ a; 
otherwise, auppcr <~ a. 

• Until auppcr ~ aiowcr < Sa, where the small positive 
constant Sa is chosen to ensure sufficient accuracy. 



> Compute the MI relay beamforming matrix as Ami = 
U*BmiU^ for each n-th subgroup, n = 1, . . . ,N. 

• Select the number of subgroup N and its corresponding 
Ami which results in the largest achievable sum-rate. 

IV. Simulation results and discussions 

In this section, the sum rate performance of the proposed 
TWR beamformers are presented. For simplicity, the channel 
correlations between the j-th and the fc-th sources are set to 
be equal, that is, Ih^^h^p = p, =^ k. 

In Fig. [T] two pairs of single-antenna source {K = 4) 
and a multi-antenna relay {M = 8) are considered in the 
simulations. The transmit power at each source Si, S2, S3 
and 54 are fixed as pi = P2 = Ps = P4 = lOdB, and the 
relay transmit power is set to be pR = lOdB. Both sources 
within each pair are set to have identical SINR requirement. 
With no channel correlation between different sources, the 
proposed MI beamformer achieves the same achievable rate 
region as the optimal MP beamforming scheme. However, 
when the channel correlations between different source pairs, 
and within each source pair, are set to be higher, the proposed 
MI beamformer achieves a smaller rate region (not shown 
here due to space constraint). This is because when the 
channels of different source pairs are correlated, the desired 
signal components that point in the direction of the inter-pair 
interference will also be suppressed by the MI beamformer, 
which in turn reduces the achievable rate for each source. 
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Fig. 1 . Compai'ison on the achievable rate region for the proposed MP and 
the MI relay beamformers with M = 8, K = 4, fixed pi = P2 = P3 = 
P4 = 10, Pji = 10, and p = 0. 

In Fig. m we present the achievable sum-rate versus SNR 
at the relay, which is defined as pn/a^, for the proposed MI 
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Fig. 2. Comparison on the sum-rate versus SNR for the proposed MI relay 
beamformer with M = 8, K = 8, p = 0. 



beamformer with various number of subgroups. Four pairs 
of single-antenna source (K = 8) and a multi-antenna relay 
(M = 8) are considered, and the channels between different 
sources are uncorrelated. Time division is used to serve dif- 
ferent subgroup. We observe that either the smaller subgroups 
or large group transmission does not always perform the best 
in TWR. With low SNR at the relay, the case of four groups 
(each with one source pair) performs the best, whereas for 
large SNR, the case of a single group (with four source 
pairs) performs the best. This is because with either large 
number of source pairs or small relay transmit power, the 
MI beamformer tends to suppress the interference more for 
the case of single group (with more interfering users), which 
results in SINR loss. 

To overcome this shortcoming, joint grouping and 
beamforming scheme can be used to reduce SINR loss as 
follows. For small SNR, we should apply the proposed 
beamformers to four different subgroups, and for large SNR, 
we should apply the proposed beamformers to a single 
group (with four source pairs). The improved sum-rate 
performance by joint grouping and MI beamforming scheme 
is highlighted by the solid line in Fig. ID The optimal 
grouping and selection of the SNR thresholds correspond to 
different number of subgroups are interesting subjects for 
further study. 

V. Conclusions 

New optimal and suboptimal beamformers for ANC-based 
TWR with multiple source pairs are derived by taking into 
account the tradeoff between the desired signals and the 



inter-pair interference. For low SNR, a better sum-rate 
performance can be achieved by first diving a large number 
of source pairs into smaller subgroups, and then apply 
beamforming to each subgroup using time division. 
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